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We derive a Bell-type inequality for observables with arbitrary spectra. For the case of continuous 
variable systems we propose a possible experimental violation of this inequality, by using squeezed 
light and homodyne detection together with methods of quantum-state reconstruction. It is shown 
that the violation is also possible for realistic detection efficiencies. 
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The theory of Bell inequalities has a rather long his- 
tory. The theory itself is named after Bell's work [1J, 
where he proposed a quantitative approach to the old 
problem of the incompleteness of Quantum Mechanics, 
posed in the famous EPR (Einstein, Podolsky, Rosen) 
paper [2]. The first result of what is now called Bell- type 
inequality is the CHSH (Clauser, Home, Shimony, Holt) 
inequality [5]- For the first time experimental violations 
of Bell inequalities were observed in [4] . Since that time 
there appeared a lot of other results, but the progress 
was mainly restricted to finite-dimensional systems and 
observables with discrete spectra. 

The literature on Bell inequalities is quite extensive, 
both the theoretical and experimental one, so here we re- 
strict ourselves mainly to some important results which 
are relevant in the context of our approach. The CHSH 
inequality was generalized to the so-called Mermin mul- 
tiqubit inequalities [5]. The Mermin inequalities where 
further extended to a complete set of inequalities for an 
arbitrary number of qubits [6 . Mermin inequalities are 
distinguished for their quantum violation, which is max- 
imal in this complete set of inequalities and it exponen- 
tially depends on the number of qubits. Among the re- 
cent trends, Bell inequalities for graph states 7 and ap- 
plications of Bell states to secure quantum key distribu- 
tion [5] have been studied. 

It is well known that Bell inequalities are fulfilled by 
local-hidden variable theories and their violation is a 
clear quantum effect. Recently an incompatibility theo- 
rem has been introduced [9], which considers nonlocal- 
hidden variable theories. It has been experimentally 
demonstrated that Quantum Mechanics violates non- 
local realistic theories [TU], for further generalizations of 
the incompatibility theorem, cf. 

As already mentioned, most of the studies dealt with 
observables with discrete spectra, usually dichotomic 
ones. For applications of the CHSH inequality for di- 
chotomic variables to continuous- variable quantum states 
see, for example, [T2J[T3]. Only recently some results for 
observables with arbitrary spectra have been obtained. 
The first multipartite Bell-type inequality for observables 
with arbitrary spectra has been given in [14 . In the bi- 



partite case it reads as 
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where arbitrary Hermitian operators A Xl A 2 and B x , B 2 
act on different sites. For two-qubit systems and di- 
chotomic observables, the right-hand side of this inequal- 
ity reduces simply to 4, so that it reads as 



{A X B X - A 2 B 2 ) 2 + (A X B 2 + A 2 B X ) 2 sC 4. 



(2) 



It has been proved by Uffink [T5], that any state (irre- 
spective if it is separable or not) of a two-qubit system 
satisfies this inequality. Hence the Bell-type inequality 
([lj can never be violated in this case. Moreover, for 
continuous- variable systems this inequality cannot be vi- 
olated in the case of quadrature measurements |16j . Vi- 
olations of the multipartite version of this inequality for 
quadratures have been obtained in the original work [14] . 
for the number of parties being at least ten. This shows 
that violations of the inequality (JTJ and its multipartite 
versions are not easy to observe. 

In the present contribution we introduce another in- 
equality which can be easily violated both by the simplest 
two-quibit Bell state and by continuous-variable states. 
As an example, we simulate an experiment with a two- 
mode squeezed vacuum. The violation can be demon- 
strated for realistic values of the squeezing parameter and 
of the detection efficiencies. 

Our main result is the following statement: any bipar- 
tite separable quantum state satisfies the inequality 



((A X B X - A 2 B 2 f + (A X B 2 + A 2 B X ) 2 ) 

> {A X B X ) 2 + (A X B 2 ) 2 + (A 2 B X ) 2 + (A 2 B 2 ) 2 . 



(3) 



At a first look it resembles the inequality (JT|). The dif- 
ference is that the squaring, averaging and the inequality 
sign are exchanged. The proof is based on the simple fact 
that any numbers a, b, c and d with ad — be satisfy the 
equality 



{a-d) 2 + {b + cf 



(4) 



Let us take Hermitian operators A\, A 2 , acting on one 
mode, and B\, B 2 , acting on the other one, and set a = 



2 



{A X B X ), b = {A X B 2 ), c = (A 2 Bx), d = (A 2 B 2 ). Then for 
a factorizable state, g — g A <g) g B , we have ad = be and 
the identity Q gives us the relation 

(Ai-Bi - i 2 B 2 ) 2 + (i x B 2 + A 2 B X ) 2 



= (A.B,) 2 + (A^) 2 + (A^) 2 + (A 2 B 2 ) 



(5) 



From the non-negativity of the variance of Hermitian op- 
erators we get 



{{A X B X - A 2 B 2 ) 2 ) + {{A X B 2 + A^) 2 ) 

^ (A 1 B 1 - A 2 B 2 ) 2 + (A 1 B 2 + A 2 B X ) 



(6) 



which, by inser ting Eq. {5}, proves the correctness of the 
inequality ^ for factorizable states. 

To extend this inequality to arbitrary separable states, 
we use the Cauchy-Schwarz inequality (J2k Pk x k) 2 ^ 
SfcPfe x L where {pk} is a probability distribution and 
Xk are arbitrary real numbers. This inequality expresses 
the non-negativity of the variance of a random variable. 
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Let us consider a separable state g — J2kPk(>k ® @k 
J2kPkQk and estimate the right hand sideboth of the in- 
equality We have 



{A 1 B 1 )l+{A 1 B 2 ) 



Pfe((ii-Bi 
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l + iA^l + iA.B,) 2 
„ + ... + (i 2 s 2 )?J 

i 2 B 2 ) 2 + (iiB 2 + i 2 S 1 ) 2 ) 



(7) 
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((ixSj - A 2 B 2 ) 2 + (i x £ 2 + A 2 5i) 2 ),. 



The first step is due to the Cauchy-Schwarz inequality 
and the second one is due to the fact that the inequality 
([3]) is valid for factorizable states, as it has already been 
proved. The last step simply expresses the linearity of 
the mean value, which completes the proof. 

If we set a = ai&i, b — a\b 2 , c = a 2 b\, d = a 2 b 2 , 
then ad — be and the identity Q simply expresses the 
multiplicativity of the norm of complex numbers 



\a\ + ih\ 2 \a 2 



ib 2 



|(a 1 + i6 1 )(a 2 + i6 2 )| 5 



We can write this equality for any number of factors and 
get corresponding multipartite Bell-type inequalities. We 
can go further and get Bell-type inequalities with more 
observables. Norms of the algebras of quaternions and 
octonions are multiplicative, and this property can be 
used to obtain the square identities generalizing the one 
given by Eq. Q , which can be used to derive multipartite 
Bell-type inequalities with four and eight observables per 
site. Details of this approach have been given in [17] . 
In this contribution we will concentrate on multipartite 
inequalities with two observables per site. 

After simple algebraic manipulations the inequality ([3| 
can be rewritten in the form 



A 1 B 1 



A X B 2 



A 2 Bt 



a\ 2§2 > \{[A U A 2 ][B X ,B 2 ] 



(9) 



where a 2 ^ = (A 2 ) — (A) 2 is the square of the variance of 

the operator A. The approach of [14] to derive the in- 
equality ([T]) is to ignore local commutators, which would 
give a trivial result in our case. We see that separabil- 
ity puts a more strict condition on the sum of the vari- 
ances then the one expressed by the inequality ( 10 ). The 



product of two commutators on the right hand side is a 
product of two local observables since it can be repre- 
sented as [Ai,A 2 ][Bi,B 2 ] = -i[Ai,A 2 }i[B x ,B 2 }. We will 
show that the inequality ^ can be easily violated. For 
the strength of violation, V, we use the ratio of the right 
hand side (containing commutators) and the left hand 
side (the sum of squares of dispersions) . When this ratio 
exceeds one, V > 1, then the inequality ([9| is violated. 
In such a case the maximal violation for a given state is 
the maximum of this ratio for all possible choices of the 
operators Ak and Bk- The inequality (10) shows that 
the maximal violation cannot exceed 2. Below we will 
see that this limit can be easily achieved. 

It is interesting to note that an arbitrary state of a two 
qubit system satisfies the following inequality: 



2 , 2 

A-i Bi ~T °A,B 2 



a AoBi 



1 



°\ 2 b 2 > ^\([ A ^M[Bi,B 2 \ 



(10) 

It differs from the inequality (JoJ) only by a constant factor 
of 1/2 on the right-hand side. The inequality (10 1 can 
be obtained from ^ using methods of [15] . To illustrate 
violations of the inequality (|9| , consider a two-qubit sys- 
tem. Let us take the operators Ak = a A k , Bk = cr£, 
where the normalized vectors and represent the di- 
rection along which the spin projections are measured. 
Then the inequality (]9| becomes 



I ~A ~B\1 



I* A ~B\ 2 



l~A ~B\1 



/ ~A ~B\1 



> 4 l(0'r 1 xr 2 CT Sl xs 2 > 



(11) 



For the Bell state 1$) 



= (l/\/2)(|01) + |10>) we have 
- r z s z — (f,s), where we set f = 



(8) ( r x, r y , — r z ); f is the reflection of r with respect to the 



xy-plane. Note that ri x r 2 = — fi x f 2 . For the vectors 
f 1, f 2 , Si and s 2 lying in the same plane with the straight 
angles between fi, fi and Si, s 2 the inequality (10 1 is 
satisfied with the equality sign. This implies that the 
Bell-type inequality ^ is violated by a factor of 2. 

For continuous variable states let us introduce the ana- 
logues of the Pauli operators via 

JV 



Sx,n = Ed 2fc )< 2fc + + \ 2k + !) < 2fc D' 
) 

N 

^(|2fe)(2fe + l| - |2fc + l)(2fc|), 

k=0 

2JV+1 2JV+1 

s z ,n = x; (-i) fc i*x*i. so,n= e l^> < fc l- 



s 



Y,N 



k=0 



(12) 
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FIG. 1: (Color online). Maximal violation of the inequality 
Q for the squeezed vacuum for different values of N. 



equation for the harmonic oscillator. Expressions for the 
non-normalizable solution ip m (x) are also known. In [18 
it is given as an action of a differential operator on the 
function tfo(x) and in [19) it is given explicitly in terms of 
the degenerate hypergeometric function. Here we could 
derive the explicit expression 



<p m (x) = wip m (x) er&(x) x 

V2" l m! 

k=0 ^ ' 



(16) 



The first of these operators commutes with the other 
three, which have the commutator rules \S a ,N, S/3,n] — 
2iS 7 , {a, 0, 7} = {X, Y, Z}, of the ordinary Pauli opera- 
tors, though their squares S x N = S Y n = $z n = &o,n 
are not the unity operator unless ./V = +00. In this limit- 
ing case we denote them simply as Sx, Sy and Sz- The 
operators S t ,n an d S r are defined in full analogy with 

<7 r . 

As a realistic example of a continuous variable state, 
let us consider the two-mode squeezed vacuum given by 
\^f) = S(z)\00), where the two-mode squeezing operator 
is defined as S(z) — exp(z*ab— za'b'). For this state (for 
real squeezing parameter z) we have 

(— ta,nh(2z)r x s x + tanh(2z)r y Sj, + r z s z ) 



and (S^S^jj) = l-tanh 4(Ar+1) (z). Figure [l] shows the 
numerically calculated maximal violation of the inequal- 
ity ([9]) for some finite values of N. The solid line corre- 
sponds to the case of N = +00. Note that for TV = +00 
we simply have 



= — tm\h{2z)r x s x 



tsnh{2z)r y s y + r z s z , 
(14) 

which in the limit z — > +00 coincides with (o'r^s) f° r 
the Bell state after reordering the components of r and s. 
Thus, in this case the violation tends to 2 when z — > +00. 

Let us now discuss experimental applications of our 
Bell-type inequality. Matrix elements g nm = (n\g\m) of 
a density operator g can be obtained from the quadra- 
ture distribution p v (x) = v (x\g\x) v , where \x) v is the 
eigenstate of the quadrature operator x^\x) v — x\x) v , 
according to the reconstruction formula 



1 



i(n—m)ip 



+ OG 



p v (x)f nm (x)dxdip. (15) 



The reconstruction kernels f nm (x) are given by f nm (x) — 
(ip n (x)ip m (x))' for n < m and f nm {x) = f mn (x) for 
n > to, where ip n (x) = (x\n) is the well known wave 
function of the nth Fock state in the coordinate repre- 
sentation, the normalizable solution of the Schrodinger 



Here [x\ is the floor of x, i.e. the largest integer that is 
not greater than x. 

Up to now we have discussed only single mode case, 
but in the case of several modes the reconstruction ker- 
nel is the product of the single mode kernels of the cor- 
responding modes. In the bipartite case under study, 
we have the two mode quadrature distribution defined 
via p v e{x,y) = <p{x\e(y\(>\x) v \y)e- The average value 



(Su.nS%n)> where V>,v = X, y, Z, 0, is given by 



P^faO^tp^W&pM, (17) 



'0 Jo 
where we denoted 

«+oo r-\-oo 



(13) 



p v e{x,y)F^N(x)F VfN (y) dxdy. 

(18) 

The phase dependent functions are Qx{<p) = - cos ip, 
& Y (ip) = -sini/3, &z(p) = 3>o(y) — an d the recon- 
struction kernels are defined via 



N 



Fx.n(x) = Fy N (x) = y^/ 2 fc,2fc+i(», 



F a ,N(x) 



k=0 

2N+1 
k=0 



(19) 



where u — I and uz = — 1. In an experiment one ob- 
tains data in the form of the table {xjj.fe, yij,k}, i,j — 
l,...,N p h, k — l,,..N poa , where one measures N qu 
quadrature values for each of N^ h chosen phase pairs 
(ifii,6j). The inner integrals in Eq. (17 1, denoted as 



PfiviPi 0), can be obtained by the sampling 



1 



JVC " U fc=l 

with i,j = l,...,iVph. Having these numbers we can 
calculate the outer integrals in Eq. ( 17 1 using the fast 
Fourier transform. 

Recently, a quantum noise redunction of a factor of 10 
has been achieved [50]. Using this, we can estimate the 
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the inequality. For this purpose a two-mode squeezed- 
vacuum state has been considered. The methods of re- 
constructing the needed correlation functions are pro- 
vided. 



FIG. 2: (Color online) Statistical distribution of the maximal 
violation obtained in a simulation of an experiment with the 
detection efficiencies rj — 0.9 (left) and 77 = 1 (right). 



degree of squeezing that can be obtained in a realistic 
experiment. For the squeezed vacuum for the disper- 
sion of x = Xa + xb we have ((Ax) 2 ) = e _2z . To find 
the maximal squeezing parameter z we have to solve the 
equation e~ 2z = 0.1, which gives z ~ 1.15. For example, 
for N — 5 we get a violation of ~ 1.9, which is very close 
to the maximal possible value. 

Finally, we discuss the maximal violation of the in- 
equality ([£]) with imperfect photodetector with the effi- 
ciency 77. The moments of the creation and annihilation 
operators /j, njn + P (r]) = (a n a^ n+p ) mcas measured by such a 
photodetector are /J, n ,n+p(v) — if l+p / 2 /j„. n+p , where the 
moments on the right hand side is the "true" moments, 
measured by the perfect photodetector with the ideal ef- 
ficiency 77 = 1 |21| . The density matrix elements can be 
obtained from the moments according to the following 
expression [22] : 



Qn.m{ll) 



^ (-1)' 



\/n\m\ ~t 



l^n-\-i,m-\-i (v)- 



(21) 



i=0 



The extension to the multimode case is straightforward. 
Thus, we can calculate the average values (S^ Q S^ Q } with 
the matrix elements obtained from the measurements 
with non-ideal efficiency 77. 

We have simulated an experiment by generating N qu = 
10 3 points (xij,k,Uij,k), k = 1, 10 3 with the distribu- 
tion p Vi 3 (x,y) for i,j = l,...,N pil = 61, where ip l = 
(7r/60)(i— 1), 9j = (7r/60)(j — 1). These points were used 
to calculate the quantities p^ v {(pi,9j), fi, v = X, Y,Z, 0, 
according to Eq. ( 20 1 . Then we have used the simplest 



numerical integration scheme to get the quantum mean 



values (S^ N S„ N ) according to Eq. (17 1. As an example 



we took z — 0.8 and T] = 0.9, 1.0. Fig7[2 shows the distri- 
bution of the maximal violation V for 2000 runs of the 
simulation. One can calculate the maximal violation by 
using the analytical expressions for (S^ N S^ N ) obtained 
with Eq. (21 1. Comparing with the values averaged over 



2000 runs of the experiment simulations, we found that 
in all four cases the difference is ~ 1%. 

In conclusion, we have derived a Bell-type inequality 
for arbitrary observables, which can be violated for dis- 
crete as well as for continuous-variable quantum states. 
We have simulated a realistic experimental violation of 
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